More about the Laughlin droplet 

N. Barberan 

Departament d'Estructura i Constituents de la Materia, 
Facultat de Fisica, Universitat de Barcelona, E-08028 Barcelona, Spain 

Abstract 

We compare the exact diagonalization ground state wave function (calculated without any re- 
striction) of a two-dimensional droplet in a perpendicular magnetic field with the Laughlin ansatz, 
as the number of electrons increases. The fully spin-polarized case for filling factor 1/3 without 
lateral confining potential or Zeeman effect is considered. We observe that the overlap decreases 
as the number of electrons increases. 

PACS numbers: 73.21.La, 73.43.Lp. 
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Most of the knowledge and understanding of the fractional quantum Hall effect states 
is due to the existence of Laughlin's wave function 1], a simple and intuitive analytical 
expression that matches very well with the ground states (gs) of the system characterized 



by the filling factor z/ = where q is an odd integer. T 



lis has been tested by numerical 



calculations of finite systems for different geometries as well as experimental 

results related with measurable quantities which can be derived from its properties, as is the 
case of the fractional charge of some of the excited states ^. We analyze the overlap between 
this ansatz and the gs wave function that comes out from numerical exact diagonalization 
of the Hamiltonian for finite systems. We also consider other properties of the gs. 

Laughlin's wave function for z/ = 1 /g is given by: 



^L = \{{z.-z,y e-^\^^\''' (1) 

where z = x + iy (m. units of = \J'^^ ^ being the magnetic field) and where the 
symmetric gauge has been assumed. The physical meaning is easily deduced from its form: 
the minimization of the energy in the partially filled lowest Landau level is obtained when 
each electron sees zeros of order q at the positions of the other electrons in such a way 
that there are non free zeros which would increase the energy. It can be verified that if 
the polynomial part is expanded, one obtains a combination of Slater determinants built up 



from Fock Darwin 
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single particle wave functions of the type: 



■'sp 



~ Z 



m g-|.|2/4 (2) 



(solutions of the non-interacting system) where m is the single particle angular momentum. 
The expansion is a homogeneous polynomial of degree M (the total angular momentum) 
and contains single particle angular momenta up to a maximum given by, 

rriL = qiN - 1) (3) 

where is the number of electrons. In order to analyze the evolution of the wave function 
given by Eq.(l) with A^, we consider the cases = 2, 3 and 4 keeping the condition given 
by z/ = 1/3 . We compare Eq.(l) with the wave function that comes out from the exact 
diagonalization (we only considered the fully spin polarized case). In the exact calculation 
no Zeeman or kinetic energy contributions are considered, electrons are confined within a 
disk due to the restriction on the total angular momentum. 



For N — 2 the dimension of the subspace characterized by the well defined quantum 
numbers M — 3, the angular momentum needed to obtain u — 1/3 for N — 2 (M = 
qN{N — l)/2) and total spin along the z-direction — 1 is equal to two, that is to say, 
the exact gs is a combination of two Slater determinants in which the largest possible single 
particle angular momentum is given by, 

,^^^^^_ (^- 1X^-2) 

or niex = 3 in this case, in agreement with the value obtained from Eq.(3). Furthermore, the 
linear combination of Slaters is the same in both cases (in Eq.(l) and in the diagonalization) 
and consequently the overlap is equal to one. That is to say, for = 2 the Laughlin wave 
functiuon is exact. For N = 3 (M = 9) the situation is not the same: the number of Salters 
involved in the g.s. is not the same, the maximum m is also different and as a consequence 
the overlap is lower than one. This tendency is increased for N — A. A summary of the 
results is given in Table 1 and Figs.l and 2 below. Nl and Nex are the total number of Slater 
determinants included in the Laughlin and exact ground states respectively, i.e., N^x is the 
dimension of the Hilbert space characterizad by M and 5*2. Pl and Pex are the normalized 
weights of the most important Slater within the expansion of the g.s. respectively. This 
Slater determinant contains, for all the values of N considered, a packet of successive single 
particle angular momenta separated from the center of the dot in such a way that they 
produce a compact ring. Remarkably this compact structure made up of successive single 
particle angular momenta wavefunctions ( m = 2, 3 and 4 for = 3) differs from the 
structure intuitively suggested by that of electrons surrounded by magnetic quantum fluxes. 
Namely, the zeros at the electronic positions seem to emphasize the short range character 
of the interaction rather than the separation between electrons. Fig.l shows the overlap 
S —\ {^L I '^ex) P for = 2, 3 and 4 and Fig.2 shows the values of Pl and Pex for different 
N. The tendency seems to indicate that the overlap will decrease as N increases and so, 
as the function given by Eq.(l) is a especially good approximation for a low number of 
electrons, it must be taken with some care for large iV in a finite system. 

It must be emphasized that the tendency of the overlap to worsen as the number of 
electrons increases would be a trivial result if the gs were a single Slater determinant and 
also if we used differet single particle wavefunctions: the exact wavefunctions for '^^x and 
approximate wavefunctions for In this case, the overlap would scale as (1 — e)^ e being 
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a small number for good trial functions. However this is not our case: the gs is a linear 
combination of several Slater determinants and in addition, we use the same Slaters to 
build up the linear combinations in each case, that is to say, the overlap depends on the 
coefficients, namely, on the electron-electron interaction and is not directly related to A^. 
As a consequence, for example, Laughlin's ansatz is the exact result for filling factor u = 1 
independently of A^, since only one determinant is involved. 

Special attention must be paid to the increasing difference between Nl and Nex due to 
the fact that the ratio between the weights of the Slater determinants that are lacking in 
related to the weights of the Slaters which are included is not negliglible. As an example, for 
A^ = 4 one of the Slaters lacking in the has a normalized weight of 0.002 compared with 
that of the most important one which is 0.331. Furthermore, loses the small and the 
large values of m and so the center as well as the edge of the droplet are poorly reproduced. 
As an example for A^ = 5, there are 27 Slaters that contain m = and only 16 of them 
are included in the expansion of "^l- However, in spite of the fact that the overlap of the 
gs wavefunction and Laughlin's ansatz worsens as A^ increases, it is not the case for some 
expected values, especially for the energy. If H has only the Coulomb contribution (as is 
appropriate for a fully spin polarized system in the lowest Landau level regime), and if we 
define the discrepancy D as [{Eex — Eij/Eex] x 100 where E^x = ex \ H \ ex) and 
£'i = I I \E'l) the results obtained are: 

For N = 2, D = 0, D = 10.5 for A^ = 3 and D = 1.3 for A^ = 4, i.e. improving the result 
as A^ increases. 

There are several previous studies on the Laughlin wave function for finite systems which 
devoted special attention to the study of the edge states: Mitra and MacDonald have 
analyzed the angular momentum distribution function for a droplet and found that the 
occupations are peaked at the edge (for A^ = 15, 20 and 25) and that it has a rapid decline. 
We believe that it can be a consequence of the reduced base implied in as it loses the small 
and large values of the single particle angular momenta. Tsiper and Goldman compare 
the density of a droplet obtained from Eq.(l) and from exact diagonalization for A^ = 5 to 
12 and u = 1/3. They obtain important differences at the centre of the droplet and a nearly 
exact coincidence at the edge. The difference in the electron-electron interaction implied 
in the Laughlin ansatz (short range interaction) and in the exact calculation (Coulomb 
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interaction) is invoked in order to explain the difference and the formation of striped states. 
We beheve that their result could be a consequence of the procedure used to obtain the 
results: within the exact calculation they truncate the base of the Hilbert space, i.e., the 
m'^^ considered is obtained by increasing Hilbert space until overlap 5* converges to at least 
three significant digits. However this procedure forces a precise coincidence with at 
the edge, giving no information about the total weight of the rest of the members of the 
base of the Hilbert space, which is not necessarily negligible, as we mentioned previously. 
As an example, for = 12 they consider m^^ = 35 (m^^ = 33) while = 143 is the 
exact single-perticle maximum angular momentum involved. At the center, their numerical 
calculation contains all the Slaters without restriction and they obtain strip-like oscillations 
on the radial electron densities which has been proved to be responsible for the observed 
unexpected behavior of the current-voltage power law 

Finally we conclude that for finite systems the overlap between the Laughlin wave func- 
tion and the exact results is worse as A^ increases. This result is similar to that obtained 
previously by Haldane ^ for spherical geometry that mimics a two dimensional homoge- 
neous system. His conclusion was that for A^ = 3 the Laughlin type function is the exact 
solution ( even for Coulomb interaction) but it is not for A^ > 4. In a recent paper by 
Yannouleas and Landman [l^] a systematic study of a system of 6 electrons in a range of 
filling factors from z/ = l/5toz/=l/9is reported. They conclude that the analytical 
model of collectively rotating electron molecules (REM) jl^ provides better representation 
of the system. Other references 0, |^ have tested some results obtained with Eq.(l) for 
some particular values of A^, however our aim is to study the tendency as A^ increases. 
The contribution of our report refers to the edge as well as the central properties of finite 
systems which can be appreciably different from those properties obtained by the use of 
the Laughlin droplet for large A^ and as a consequence, the differences can be significant at 
the thermodynamic limit. However, the evolution of the overlap of the wave functions does 
not necessarily characterize the evolution of the expected values of some operators as was 
previously pointed out for the eneregy operator. 
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been performed under Grants No. BFM2002-01868 from DGESIC (Spain) and No. 2001GR- 
0064 from Generalitat de Catalunya. 
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FIG. 1: Overlap between the Laughlin wave function and the exact diagonahzation result as a 
function of N 
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FIG. 2: Evolution of the normalized weight of the most important Slater determinant within the 
wave function expansion as N increases 
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Table 1: Fully spin polarized droplet for Sz = N/2 and v = 1/3. A4 is the 
total angular momentum, m.j^ and nipx are the maximum single particle angular 
momentum for each M . and N^^ are the number of Slater determinants 
involved in the expansion of the ground state. Pj, and P^^ are the normalizad 
weights of the most important Slater within the ground state ; for TV — 2 it 
contains: m-\ = 1 and ma = 2 (let's say (L2)), Af = 3: (2,3,4), TV ^ 4: 
(3, 4, 5, 6), A' = 5: (4, 5, 6, 7,8). S is the overlap (see text). 
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